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ABSTRACT 

We present Hubble Space Telescope + Cosmic Origins Spectrograph (HST+COS) 
data of the eclipsing double white dwarf binary CSS 41177. Due to the tempera¬ 
ture difference between the two white dwarfs, the HST+COS far-ultraviolet data are 
dominated by the hot, primary white dwarf and allow us to precisely measure its tem¬ 
perature (Ti). Using eclipse observations, we also tightly constrain the temperature 
of the cooler secondary white dwarf (T 2 ). Our results, where T\ = 22439 ± 59 K 
and T 2 = 10876 ± 32 K, with the uncertainties being purely statistical, place the 
secondary inside and close to the blue edge of the empirical instability strip for low 
temperature hydrogen-atmosphere white dwarfs. Dedicated high-speed photometry is 
encouraged to probe for the presence of pulsations, which will constrain the border of 
the instability strip as well as probe a new region of low gravity within the strip. 

Key words: stars: individual: SDSS J100559.10+224932.3 - white dwarfs - binaries: 
eclipsing - techniques: spectroscopic. 


1 INTRODUCTION 

CSS 41177, also known as SDSS J100559.10+224932.3, was 
discovered to be an eclipsing binary containing a white dwarf 
during a search for transiting planets around white dwarf 
stars using Catalina Sky Survey (CSS) data dDrake et all 
l20ldl . Using an additional detection in the 2MASS survey, 
the authors put a weak constraint of spectral type M6 on 
the secondary star. However, given that the CSS light curve 
was not well sampled, they also mention that the data are 
consistent with a small faint object as the companion to the 
whi te dwarf. _ _ 

IParsons et al.l (1201 il l obtained highe r quality photome - 
try using the Liverpool Telescope (LT, ISteele et al.1 120040 
+ RISE camera, which led to the detection of the sec¬ 
ondary eclipse and the realisation that CSS 41177 is a dou¬ 
ble white dwarf binary with an orbital period of 2.78 hrs. 
CSS 41177 was the second eclipsing double white dwarf bi¬ 
nary to be discover ed, and it is still one of only five known bi- 


naries of this type (llSteinfadt et al. 

20ld;IVennes et al] 20111; 

Brown et al.l 201ll: Kilic et al. 2011 

). In addition, CSS 41177 


is the only one of these that is a double-lined spectroscopic 
binary, making it ide ally s uited for a detailed parameter 
study. IParsons et al.1 (1201 il l combined the phase-resolved 
LT + RISE photometry with phase-resolved Gemini Multi- 
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Object Spectrograph (CMOS) spectroscopy to determine or¬ 
bital parameters and the white dwarf masses and radii. 


In January 2012, higher signal-to-noise photometric 
data were obtained using the high-speed camera ULTRA¬ 
CAM ([Dliillon et al.l l2007 1. which images in thr ee bands si ¬ 
mult aneously (Sloan Digital Sky Survey (SDSS lYork et al.1 
I2000I1 u'g'r' filters in this case), while it was mounted on 
the William Herschel Telescope and New Technology Tele¬ 
scope. Combined with Very Large Telescope (VLT) + X- 
Shooter spectroscopy with a time resolution of 5-6 minutes, 
the unce rtainties on the ma sses and radii were substantially 
reduced dBours et al.ll2014l l. In addition, the changes in the 
relative eclipse depths in the three ULTRACAM bands con¬ 
strained the primary white dwarf’s temperature, while the 
differences in primary and secondary eclipse depths also con¬ 
strained the secondary’s temperature: Ti = 24407(654) K 
and T 2 = 11678(313) K. Combined with the surface gravi¬ 
ties from the ULTRACAM modelling, it was found that the 
secondary, cooler white dwarf was positioned very close to 
the empirical ZZ Ceti instability strip. Most white dwarfs 
found in this region of temperature - surface gravity param¬ 
eter space show pulsations on timescales of 7-70 minutes. 
However, the ZZ Ceti strip is largely defined by a group of 


pulsating white dwarfs with surface gravities log g ~ 8 ± 0.5 
dGiannirias et al1l200f ) and a hand ful of pulsators near log g 
= 6.5 ( Hermes et al.l 20121 . 12013bl ). Because the cool, sec¬ 
ondary white dwarf in CSS 41177 lies in between these two 
groups, near log g = 7.3, confirmation of its (non-)pulsating 
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Figure 1 . HST+COS light curve of CSS 41177, displayed in 5 sec¬ 
ond bins, showing the eclipse of the hotter white dwarf (primary 
eclipse) and the model used to determine the mid-eclipse time 
(see Section ej. The data between the dotted lines was excluded 
from the spectroscopic fits. 


nature could help determine if the instability strip is continu¬ 
ous and if the current extrapola ted empirical boundaries are 
correct, as predicted by theory (IVan Grootel et aLll2013lh In 
our fits of the optical ULTRACAM light curves, the values 
of Xi and T 2 were correlated, but at far-ultraviolet wave¬ 
lengths the hotter primary white dwarf dominates the flux. 
Here we describe an improved parameter determination of 
CSS 41177 based on far-ultraviolet HST+COS observations. 


2 HST+COS OBSERVATIONS 

We obtained far-UV HST+COS dGreen et al.ll2012h data of 
CSS 41177 as part of Cycle 21, program ID 13421. CSS 41177 
was observed during two consecutive HST orbits on May 28, 
2014. Each orbit was split in two exposures to retain a good 
signal-to-noise ratio while minimising the effects of fixed 
pattern noise using all four FP+POS settings. The expo¬ 
sures obtained were 19, 20 and 2x24 minutes long and each 
used the G140L grating at a central wavelength of 1105 A. 
During the third exposure a prima ry eclipse o c curre d 
(cycle 7491, using the ephemeris from [Bours et all J2014lb 
with a mid-eclipse time Imjd(utc) = 56805.41665(34), 
equivalent to the barycentrically-corrected time 
^bmjd(tdb) = 56805.41628(34), see Fig. 0 

We removed the part of the exposure in which the 
eclipse occurred for the following analysis using calcos to 
filter out the data at 900 - 1100 s. We also removed geo¬ 
coronal Lyman-a and Ol (1304 A) emission and interstellar 
absorption lines (see Section [4]). The lower limit for our data 
is set at 1150 A, to exclude artificial features near the edge of 
the observable wavelength range. The upper limit is fixed at 
1700 A, motivated by the decreasing sensitivity of COS and 
the therefore increasing difficulty of relative flux calibration 
at longer wavelengths. 


3 RESULTS 

3.1 Finding temperatures and surface gravities 

To determine the whi te dwarf temp eratures, we fit the data 
with model spectra dKoesteij l2010i b with a mixing length 
ML2/a = 0.8. However, simply fitting the HST+COS data 
with a single white dwarf model spectrum over predicts the 
flux from the hot white dwarf at optical wavelengths. Be¬ 
cause the secondary white dwarf starts contributing to the 
total flux at A > 1400 A, the slope of the HST+COS spec¬ 
trum is shallower than would be the case if only the primary 
white dwarf is visible in this regime. As a result, the temper¬ 
ature we would find for the primary white dwarf is too low. 
Therefore we use an iterative process which, after the first 
step, includes the secondary white dwarf’s contribution to 
the far-UV flux in order to find a more accurate temperature 
for the primary white dwarf. 

In the first step we fit the HST+COS spectra with a 
single white dwarf model spectrum using a Markov-chain 
Monte Carlo (MCMC) analysis usin g the python package 
emcee dForeman-Mackev et al.ll2013l ). This MCMC fit in¬ 
cludes as free parameters the primary white dwarf’s tem¬ 
perature (Ti) and surface gravity (log gi), the reddening 
towards CSS 41177 ( E(B — V)), and a scale factor to allow 
for the unknown distance to the binary. The surface gravity 
is constrained by a Gaussian prior at log g\ = 7.32 ± 0.02 
ll Bours et al .11201 il ) and the red dening by a Gaussian prior a t 
E(B - V) = 0.0292 ± 0.0009 llSchlafly fc Finkheinerll201lfl . 

We then fitted the L1LTRACAM li ght cu rve data as pre¬ 
sented and discussed in iBours et alj (l2014h . but now with 
a prior on the temperature of the primary white dwarf 
Ti based on the results from the MCMC analysis of the 
HST+COS data described above. Briefly, the ULTRACAM 
data include twelve primary and nine secondary eclipses, 
all observed simultaneously in the SDSS u'gV filters. The 
eclipse depths in these three bands constrain the relative 
temperatures of the white dwarfs. The free parameters in 
this fit are T\ and X 2 , the radii relative to the binary’s semi¬ 
major axis Ri/a and i? 2 /a, the radial velocity amplitudes 
K 1 and K 2 (constrain ed by X-Shooter spectroscopy also pre¬ 
sented in lBours et al.ll2014f) . the zero-point of the ephemeris 
To, the orbital period and inclination P or b and i, the redden¬ 
ing E(B — V) and an offset S from phase 0.5 of the secondary 
eclipse that allows for a Rpmer delay or eccentricity of the 
orbit. The tight prior on Ti given by the HST+COS data 
leads to a much more precise temperature for the secondary 
white dwarf X 2 as well. 

This, in turn, allows us to refit the HST+COS data 
while accounting for ‘contamination’ to the UV-flux by the 
secondary white dwarf by subtracting a reddened model 
spectrum from the data, with values T 2 and log g 2 fixed at 
the mean values from the ULTRACAM light curve fit. The 
new value found for T\ then becomes a more accurate prior 
in the light curve fits. By including the contribution of the 
secondary white dwarf, the primary white dwarf’s tempera¬ 
ture increases by ~ 150 K after the first iteration, and ~ 10 
K after the second iteration. The latter is well within the 
statistical uncertainty and eliminates the need for further 
iterations. 

The converged results are listed in Table 0 We have 
experimented with placing the long wavelength cut on the 
HST+COS spectra at 1500 A, which decreases the primary 
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Figure 2. HST+COS spectrum of CSS 41177 (grey dots, 1 A bins; crosses denote data excluded from fits). The solid black line is 
the sum of the primary and secondary white dwarf model spectra corresponding to the best fit from the HST+ULTRACAM analysis 
described in Section l3.ll Also shown is the GALEX far-UV flux (yellow dot, errorbar too small to be seen) and the position of interstellar 
absorption lines (vertical grey dotted lines). 



wavelength (A ) 


Figure 3. Top panel: model spectra for the double white dwarf binary CSS 41177 and individual white dwarfs (dotted grey lines). The 
HST+COS and SDSS spectra are shown in dark grey. Solid yellow dots indicate the GALEX far-UV, GALEX near-UV and SDSS ugriz 
fluxes (errorbars too small to be seen). Bottom panel: residuals of the model CSS 41177 spectra folded through the SDSS filter curves 
with respect to the measured SDSS fluxes. 


white dwarf’s temperature by ~ 150 K, roughly illustrating 
the magnitude of the systematic uncertainties. The uncer¬ 
tainties in Table Q] should be considered as statistical uncer¬ 
tainties only. 

The HST+COS data and the best synthetic model for 
the double white dwarf from the MCMC analysis are shown 
in Fig. [2] Fig. [3] shows the model spectra for the individ¬ 
ual white dwarfs as well, and covers the far-UV through to 
the near-IR wavelength range. The CSS 41177 models are 
the sum of two white dwarf model spectra, accounting for 


the difference in their surface area. The parameters (T e fr, 
log g) for the primary white dwarf’s model spectrum are 
those from the best fit in the last iteration to the HST+COS 
data, and the parameters for the secondary white dwarf are 
taken from the corresponding model to the ULTRACAM 
light curves. 

The CSS 41177 model spectrum is in good agreement 
with the HST+COS data as well as with the SDSS data at all 
wavelengths. The slight over- and under-predictions of the 
(AB corrected) flux that can be seen in the residuals in the 
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Table 1. White dwarf parameter results from the MCMC anal¬ 
yses performed on the HST+COS and ULTRACAM data. For 
comparison, the results from [Bours et al.l d2014f) based on the 
ULTRACAM light curves alon e are shown as well ( E[B — V) was 
fixed to the mean value from ISchlegel et al .111 998 *1. Numbers in 
parentheses indicate uncertainties in the last digit(s). 


parameter 

HST+ULTRACAM 

ULTRACAM 

Ti (K) 

22439(59) 

24407(654) 

T 2 (K) 

10876(32) 

11678(313) 

E(B - V) (mag) 

0.0292(9) 

0.0339 (fixed) 

log(ffl) 

7.322(15) 

7.321(15) 

l°g(S2) 

7.305(11) 

7.307(11) 

Ri/a 

0.02508(22) 

0.02510(22) 

R 2 /a 

0.02340(28) 

0.02332(28) 

minimum x 2 

44465 

44457 

degrees of freedom 

37745 

43351 


bottom panel of Fig. [3] amounts to 2-3 times the uncertainty 
on the measured flux or, equivalently, a few %. Because the 
SDSS uncertainties do not include systematic uncertainties, 
we do not think this difference is significant. 



Figure 4. ZZ Ceti diagram, with the filled, black star showing the 
position of the secondary white dwarf in CSS 41177 at the overes¬ 
timated values that would have been recovered if the parameters 
were found by fitting the seco ndar y white d warf’s spectra with ID 
model atmospheres dTremblay et a l. l2013h . The open star indi¬ 
cates the direct values found from our fit — see text. The black and 
grey dots indicate confirmed pulsating and non-pulsating white 
dwarfs respectively, and the dotted lines are em pirical boundaries 
for the instability strip llGianninas et allllOl l l. 


4 METAL LINES IN THE HST+COS SPECTRA 

The HST+COS spectra reveal a handful of absorption lines. 
The position of these lines coincides with known interstel¬ 
lar absorption features: Sill 1260 A; Ol 1302,1304 A; Cll 
1334/5 A. Given that the lines also do not shift position be¬ 
tween the four individual exposures we conclude that they 
are indeed of interstellar origin. The Sin 1265 A excited state 
is not detected, further corroborating the interstellar nature 
of the metal lines. 


5 A PULSATING SECONDARY WHITE 
DWARF? 


Depending on their atmospheric composition, white dwarfs 
experience non-radial gravity-mode pulsations at different 
epochs in their cooling curves. Based on extrapolating obser¬ 
vations of white dwarfs with lower and higher surface grav¬ 
ities, this instability should occur close to Twd = 11000 
I< for white dwarfs with a mass of Mwd ~ 0.3 and 
hydr o gen-rich atmospheres (DA) (see e.g. iBergeron et al.l 
1 19951 : IGianninas et al .1 1201 il l . In the log g - T e ft plane this 


region is called the ZZ Ceti instability strip. 

The values of the secondary white dwarf’s temper¬ 
ature and surface gravity from both o ur analyses are 
shown in Fig. HI where known p ulsato rs llGianninas et aD 
l201ll : IHermes et al.ll201^ . l2013al la. [201^ P^rzas_et - ^lJ 20_l_o|) 

l201ll : 


and non-pulsators a t 10 mmag llGianninas et alJ 20 111 : 
ISteinfadt et alj 120121 : IHermes et al.l I2012L l2013bl cj)" are also 
indicated. The values used for low-mass white dwarfs are 
the up dated atmospheric parameters from IGianninas et al.l 
d20141 ). who used the latest model spectra, with a 1- 
dimensional mixing length theory to approximate convec¬ 
tion, to fit the white dwarf Balmer lines. It has been known 
for some time that fitting Balmer lines with ID models over 
predicts surface gravities and temperatures in the regime 


where w hite dwarf atmospheres beco me convective (T e g < 
13000 K: iTremblav et al.ll201lll20Hl) . 

The method presented in this paper does not rely on fit¬ 
ting Balmer lines, but rather uses the continuum of the spec¬ 
tral energy distribution and the L1LTRACAM light curves 
to find the temperature and surface gravity of the secondary 
white dwarf in CSS 41177. These results are likely close to 
the true atmospheric parameters, and therefore presumably 
similar to results that would be obtained if it was possible 
to fit the secondary’s Balmer lines directly with the latest 
3D model spectra. For this reason, we corre ct our results in 
Fig. f using the 3D to ID correction given bv lTremblav et all 
d2013l) . to facilitate direct comparison with other data points 
in this figure. The filled star indicates the position in the ZZ 
Ceti strip that would have been recovered if fits to the sec¬ 
ondary star’s Balmer lines were performed with ID model 
atmospheres. The open star shows the atmospheric param¬ 
eters found with our method, also listed in Table |T| Note 
that the average of the 3D models is performed over constant 
Rosseland optical depth, which produces very similar results 
to full 3D spectral synthesis at a ll wavelengths for thes e tem¬ 
peratures and surface gravities dTremblav et al.ll20111 ). 

Our results place the secondary white dwarf inside the 
empirical instability strip, just under 350 K from the blue 
edge. With the available ULTRACAM data, pulsations with 
a relative amplitude exceeding 0.5% in the light of the sec¬ 
ondary (i.e. correcting for the 80% contribut ion of the hot 
primary) were ruled out dBours et al.ll20L4h . However, as 
noted in that paper, the ULTRACAM data were not ideally 
suited for finding pulsations. The observations targeted the 
eclipses, and mostly covered only brief stretches (5-10 min¬ 
utes) of out-of-eclipse dat a. Typically, pulsation perio ds are 
between ~ 7-70 minutes (IHermes et al.l I20T2I . l2013bl ) , and 
increase with decreasing surface gravity and temperature. 
Therefore the ULTRACAM observations are too short. Fur- 
therm ore, pulsating white d warfs can have amplitudes < 
0.5% (iMukadam et al.ll2006l f. This might explain why we 
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have not been able to detect any pulsations. New, dedi¬ 
cated observations will reveal whether pulsations are indeed 
present or not. 

Note that we have assumed a prior for the reddening at 
E(B—V) = 0.0292 ± 0.0009 in our MCMC analyses, which is 
the total estimated reddening in the direction of CSS 41177. 
Since the reddening affects the flux at shorter wavelengths 
more than the flux at longer wavelengths, a smaller E(B — 
V) results in a shallower slope of the HST+COS spectra. 
This translates to a lower temperature for the primary white 
dwarf, and, through the ULTRACAM light curves, also in a 
lower temperature of the secondary white dwarf. Therefore 
any decrease in the reddening will push the secondary white 
dwarf further into the instability strip. 


6 CONCLUSIONS 

We have presented HST+COS spectra of the eclipsing dou¬ 
ble white dwarf binary CSS 41177. By fitting the far-UV 
data with model spectra, these allow us to precisely mea¬ 
sure the primary white dwarf’s temperature: T\ = 22439 ± 
59 ± 150 K. Combined with the eclipse data in previously 
presented u'gV ULTRACAM data this puts a very tight 
constraint on the secondary temperature as well: T 2 = 10876 
± 32 ± 150 K. 

The ULTRACAM light curves, combined with con¬ 
straints fro m X-Shoote r tim e-series spectroscopy also 
presented in lBours et al.l (l2014l l. also constrain the surface 
gravity: log g 2 = 7.305 ± 0.011. These results place the 
secondary white dwarf inside the pulsational instability 
strip for white dwarfs with hydrogen-rich atmospheres, and 
in particular in the gap between the canonical log g ~ 8 
white dwarfs and the extremely low mass white dwarfs 
at log g ~ 6.5. The previously presented ULTRACAM 
data include only short stretches of out-of-eclipse data, 
in which we were unable to detect pulsations. Given the 
results presented here an d the predictions made by theory 
llVan Grootel et all 1201. 'll ), dedicated high-speed photomet¬ 
ric observations resulting in higher signal-to-noise ratio 
data may still reveal the presence of pulsations. This white 
dwarf is particularly interesting because it is positioned in 
an as-yet-unexplored part of the ZZ Ceti diagram, as well as 
being close to the edge of the instability strip. Determining 
if it pulsates or not will help determine where the border 
of the instability strip is and whether the strip is con¬ 
tinuous from high mass to extremely low mass white dwarfs. 
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